/?-decay properties for neutron-rich Kr-Tc isotopes from deformed pn-QRPA 

calculations with realistic forces 
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In this work we studied /3-decay properties for deformed neutron-rich nuclei in the region Z=36-43. 
We use the deformed pn-QRPA methods with the realistic CD-Bonn forces, and include both the 
Gamow- Teller and first-forbidden types of decays in the calculation. The obtained /3-decay half-lives 
and neutron-emission probabilities of deformed isotopes are compared with experiment as well as 
with previous calculations. The advantages and disadvantages of the method are discussed. 
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I. INTRODUCTION 

Decay properties such as the half-lives and /3-delayed 
neutron emission probabilities are important inputs for 
the simulations of the r-process nucleosynthesis which is 
believed to be responsible for the production of heavy 
elements in our universe. In order to understand the 
observed elements abundance, one needs measurements 
together with models for making accurate predictions for 
these global properties of atomic nuclei out to the neu- 
tron drip-lines, especially those neutron-rich nuclei along 
r-process paths [if, Q • 

Recently, a group from RIKEN has performed a se- 
ries of half-life measurements for the neutron-rich Kr-Tc 
isotopes 0. For some of these nuclei, some differences 
from the previous measurement has been found, while 
for others, the half-lives are measured for the first time. 
These measurements give us more information for exotic 
neutron-rich nuclei and also offer us more information 
relevant for the r-process flow path around the A=130 
peak. 

These new measurements serve as good tests or con- 
straints for theories. The theory for such calcula- 
tions from gross to microscopic have been developed for 
decades. There have been global estimations of half-lives 
from the gross theories such as those in Ref. [5] which 
treats the half-lives as functions of the Q-values, pro- 
ton (Z) and neutron (N) numbers. More microscopic 
methods have been developed in Ref. [l| with deformed 
pn-QRPA methods for the Gamow- Teller (GT) type de- 
cays with the residual interaction from the phenomeno- 
logical pn forces in J v = 1 + channel. Because of the 
phenomenological nature of the forces designed only for 
the GT channel, the First-Forbidden (FF) part was esti- 
mated by gross theory instead of direct calculations [l[ . 
Despite these approximations, Ref. |I| gives a reason- 
able average agreement with half-lives over the nuclear 
chart and for the new RIKEN data, but there are order 
of magnitude deviations between experiment and theory 



m many cases. 

As an improvement, we proposed the adaption of re- 
alistic forces for the residual interactions. This method 
was first developed by the Tubingen group H, [l(| for 
/3/3-decay. The advantage of this method is that we are 
not restricted to GT channels. With just two parameters 
for the residual interactions, we can calculate all possible 
decay channels such as those transitions from + ground 
states to 0~,1~,2~ final states. This gives explicit re- 
sults from microscopic calculations for FF decays which 
in some nuclei may play an important role. On the other 
hand, inclusion of all spin-parities can give us better de- 
terminations of the spectra of odd-odd nuclei that makes 
the decay energy required for phase-space factors much 
more accurate. Both are advantages which can make a 
better prediction for /3-decay properties. 

This article is arranged as follows; in Sec. II we in- 
troduce the formulations of the calculations for /3-decay 
properties and the many-body approach we adopt, in Sec. 
Ill the results and comparison to experiment and other 
models are discussed. The Conclusions are given in Sec. 
IV. 



II. FORMALISM 



The half-lives of /3-decaying isotopes can be expressed 



as: 
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where T is the decay width and has the form: 
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Here the sum runs over all the possible states j for 
the final nuclei, fj is the phase-space factor for state 
j, which is a function of the nucleus radius R, nu- 
clear charge Z and the energies of the emitted electron 
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eo = E e /m e = (Q — Ej)/m e in unit of the electron mass 
m e , while Mj is the nuclear matrix element. In this work, 
we consider two kinds of decays, the Gamow- Teller (GT) 
and first-forbidden (FF) decays. The phase-space factors 
for both decays are expressed analytically in Ref. [4j and 
references therein. The nuclear matrix elements describe 
the nuclear-structure part of the /3-decay. They can be 
expressed in the form: 

Mj = (i|0 J |i>, (3) 

where |i) is the ground state for the parent nucleus and 
O 1 are the transition operators for /3-decay and have dif- 
ferent forms and selection rules for different types of de- 
cay /. From the expression in ^ for the half-lives for 
decays to each final state, one needs the information for 
the excitation energies of the final states and the ma- 
trix element between the initial and final states. For 
different type of nuclei, we will use different many-body 
treatments as explained below. 

Various methods have been applied to this calculation 
each with some limitations. The large-basis shell model 
can account for many of the correlations, but the num- 
ber of orbitals that can be considered is restricted by the 
computational limitations to dimensions of on the order 
of 10 10 , and in practice is used only for the Gamow- Teller 
decay in light nuclei (A < 60). All of the nuclei consid- 
ered here are outside of this range. 

Thus one needs other methods that involve various ap- 
proximations. In our case, we adopt the deformed ver- 
sion of pn-QRPA with realistic forces as first introduced 
in Ref. [9(. With the adiabatic Bohr-Mottelson approxi- 
mation, one can prove the equivalence of the calculations 
performed in the laboratory systems and intrinsic sys- 
tems Q (In the intrinsic system, the z axis is attached 
to the symmetric axis of the nucleus). Thus we will per- 
form our calculation in the intrinsic system without the 
consideration of rotations of the nucleus, and adopt the 
axially symmetric wave functions in the intrinsic frame. 
Details of the wave function calculation is described in 
Ref. [9( as well as the treatment of BCS pairing in the 
deformed nuclei. With the BCS pairing, one defines the 
Boglyubov quasi-particle creation and annihilation oper- 
ators: 




where the annihilation operators annihilate the BCS vac- 
uum, a T \BCS) = 0. Here c and are single-particle 
annihilation and creation operators, and u's and w's are 
BCS coefficients from the solutions of BCS equation. 

The intrinsic excited states \K™ , to) are then generated 
by the QRPA creation operators acting on the ground 
states [ig |: 

\K\m) = Qk,JO+,.> 

Qif\m = ^2 X^ x* A pn K „ — Y™ jK -n Apn,K^ ■ (5) 
p . 1 1 



Here the two quasi-particle creation and annihilation op- 
erators are defined as: A pnK7T = <xj,a:~ and A 
apa n , with the selection rule: K = Q, p — il n . In order 
to obtain the forward and backward amplitudes X's and 
Y's, one needs to solve the QRPA equations: 

( A{K«) B{K*) \ ( X" K \ \_ ( X" K \ \ 

\ B(K*) A(K*) J \ -Ygt J ~ K *' m \ -Y™ J ■ [0> 

The expression of matrices A and B and the details of 
the interactions are discussed in Refs. 

HQjI. The QRPA 

equations can be solved by diagonaliztion following the 
method in Ref. [12j |. The solutions contain the informa- 
tion of the energies from the eigenvalues and the struc- 
ture information from the forward and backward ampli- 
tudes. With the states constructed from the solutions of 
QRPA equations, we can calculate the beta-decay half- 
lives. Here we briefly introduce the details of the calcu- 
lations for different types of nuclei. 

First, for the decays of even-even nuclei [the first even 
(odd) refers to the proton number Z and second even 
(odd) refers to the neutron number N], the parent nu- 
clei have ground states with all the neutrons and protons 
paired. Thus it is the BCS vacuum with J v = + . The 
excited states for daughter nuclei in pn-QRPA formal- 
ism are just those we constructed above, we choose the 
states with the lowest eigenvalues to be the ground states 
of the corresponding odd-odd nuclei which are the decay 
products. (This requires a calculation over all the pos- 
sible spin projections and parities K n ). The excitation 
energies for each states are E^.m — i^>K v ,m ~ u g-s.- The 
matrix elements of the decay to the mth states with spin- 
parity K w can then be derived as following: 

= ^5 K M{p\0{i\n){X^ nK ^u p v n +Y^ K ^v p u n ). (7) 

pa 

Here (p|0j^|n) is the single-particle transition matrix 
element in the deformed basis, which can be written 
as a decomposition over the reduced matrix elements 
in spherical harmonic oscillator basis \9\: (p|Oj^|n) = 
E, p , Vn F PnlnvM\ 0l \M/V2TTT. Here for GT de- 

cay the operator has the form 0^ T — (tkt + , with the 
selection rules AK = 0, ±1 and An = 1, while the ex- 
pressions for the first-forbidden beta decay are more com- 
plicated with six components with different spin-parity as 
introduced in Ref. [TJ]. We will not give the explicit ex- 
pression for these components here, but mention that of 
which two have the selection rules AK — 0, Ait = — 1, 
three have the selection rules AK — 0, ±1, Att = — 1 and 
one has AK = 0, ±1, ±2, Att = -1. 

For odd-mass nuclei, we follow the method in Ref. 
except that the Av = case (defined in Ref. 0), wo 
consider only the single-particle transitions without the 
corrections from particle-vibration couplings for simplic- 
ity. The states of the odd-mass nuclei can be described 
as one corresponding quasi-particle (proton or neutron) 
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excitation on the even-even ground states: 

\Z±l,N,i) = al l \Z,N,0+) 
\Z,N±l,i) = al ti \Z,N,Q+), (8) 

or a pn-QRPA excitation state with one spectator single 
particle (or hole): 

\Z + l,N,im) = a^W-va) = a^Q^jZ, N, 0+) 
\Z,N-l,im) = al^-ra) = a^Q^jZ, N, 0+). (9) 

Here the label i is the spectator nucleon which doesn't 
participate in the decay process. The ground states of 
these nuclei are the one quasi-particle states with the 
lowest quasi-particle energies E Ti q; here r can be proton 
or neutron. Energies for the states in (j8]) are simply 
the differences between the quasi-particle energies E = 
Ei — Eq. While for the states in ([9]), we use different 
treatments for the energies as that in Ref. 1]: E — 
Q — Qee + E. m , E rn is the actual pn-QRPA excitation 
energy in the even-even system and Q — Qee accounts 
for the difference of the Q values between the odd mass 
isotope and the corresponding even-even one. The matrix 
elements for the Av = single-particle transitions are 
then expressed as the leading order terms 0: 

(Z + 1, N, i\0^f\Z, N + 1, 0) = u p ,i«n,o 
(Z,JV-l,»|Oj[ f |Z-l,JV,0) = up li tt n ,o(p,*|0^|n,0)(10) 

and for the Av = 1 spectator case (also defined in Ref. 
0): 

(Z + 1,N, Qm\0 T M \Z,N + 1,0) 

= 6 KM(X™n,K«UpV n + Y™ J( „V p U n )(p\0 I M \n) 

(Z,N -l,0m\O T M \Z-l,N,0) 
= J2 s KM(X™ tK „u p v n + Y™, K „v p u n ){p\0{ I \n).(ll) 

pjtp ,n 

Here contributions from the orbit occupied by the spec- 
tator particle (hole) are excluded. 

Finally, we discuss the case of odd-odd nuclei following 
the treatment from Fig. 3 in Ref. 0. In this scenario, 
the collective effect is excluded, and for the ground states 
of the odd-odd nuclei, one has simply one-neutron parti- 
cle and one-proton hole acting on the even-even ground 
states. For the daughter even-even nuclei, the ground 
states are obviously the BCS vacuum, while for the ex- 
cited states, there exists two different types: the first case 
is the two-quasiparticle excitation neglecting the collec- 
tivity, the excitation energies are: E = E p j + E Pt o (lower 
left panel in Fig. 3 of Ref. 0) or E = E n>i + E n<0 (up- 
per right panel in Fig. 3 of Ref. 0); the second case is 
the pn-QRPA excitation acting on the odd-odd ground, 
the energies are E = Q — Qee + E m (upper left panel in 
Fig. 3 of Ref. 0), Qee and E m have the same meanings 
as in the odd mass case. In the first case, the transition 
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FIG. 1: Illustrations of calculated branching ratios of the first- 
forbidden decay for different isotopes. The diagonal line in- 
dicates the approximate r-process path. 

matrix elements are just those of (fT0| with one spectator 
proton hole or neutron particle. The special case is when 
the final states is the ground states with all neutrons and 
protons paired, in this case, the excitation energies are 
instead of 2Eo (lower right panel in Fig. 3 of Ref. 0). 

The same matrix elements as in (|f f [) are adopted with 
the exclusion of contributions from the orbits occupied 
by the unpaired neutron (particle) and proton (hole) of 
the odd-odd ground states. In this approximation, the 
result is independent of the odd-odd ground state spin, 
but the lack of collectivity will give an over-estimation 
on the excitation energy of the two quasi-particle final 
states. This makes the beta-decay Q values too low, and 
increases the beta lifetime. However, as discussed further 
below, the /3-decay lifetimes for odd-odd nuclei are not 
so important for the r-process path. 

Following the definition in Ref. , the overall measure 
of error in the deviation of theory from experiment is 
defined by: 

, / ^1/2, calc x 

r = log 10 (- ) 

l l/2,exp 

s,=[-i> i ) 2 ] i/2 

i—i 

£1° = f0 Sr (12) 

Here, T,}. is defined to be the total "error". This will be 
used later for quantitative estimation on the quality of 
calculations. 



III. RESULTS AND DISCUSSION 

The single-particle wave functions and energy levels 
are obtained by solving Schrodinger equation with the 
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FIG. 2: The half-life dependence on g pp for different even- 
even isotopes. The upper panel are results with quenched 
axial vector coupling constant qa = 0.75(?ao, while for the 
lower panel the bare one <?ao = 1-26. 



deformed Woods-Saxon potentials. The parameters of 
the Woods-Saxon potential are taken from Ref. [lH • For 
the choice of the model space, based on previous experi- 
ence, we start from the Qftio to one major shell above the 
Fermi surface of either proton or neutron, depending on 
which is close to zero-energy, so for most isotopes here 
with TV < 70, the model space is — 5hui. For the de- 
formation parameters, we choose either the experimental 
ones if available (using the treatment for deformation in 
Ref. [l(|) or those predicted by Ref. [15]. For the pair- 
ing interaction, we adopt the Bruckner G-matrix with 
the Charge-dependent-Bonn force. The detailed calcula- 
tion procedure and choice of parameters are described in 
Refs. f§-[llj]. The same interaction is used as the residual 
interaction for pn-QRPA phonons. 

As mentioned in the introduction, one of the new as- 
pects of our work is the introduction of a microscopic 
calculation for the first-forbidden decays. Fig. Q] shows 
the region of nuclei covered in this paper. For each nu- 
cleus we show the branching ratio for the forbidden decay. 
They vary from several percent to at most 30% in this 
region. Thus our discussions below apply mainly to the 
Gamow- Teller decay aspects of calculations. In other re- 
gions, the FF contributions are more important and an 
accurate determination of these decay widths can give a 
better accuracy, an example is given in Q for N = 126 
isotones where, for some isotones, the first forbidden de- 
cays contribute more than 80% of the decay width. Over- 
all, an explicit calculation of FF decay is required for a 
complete account of the beta-decay in the r-process path. 

Two parameters are introduced as described in Refs. 
the renormalized particle- hole (g p h) and particle- 
particle (g pp ) strengths. However, the fitting procedure 
of these parameters is a bit different from Refs. foL fioj. 



For the renormalized particle- hole strength g p h, the usual 
way is to fit the position of the Gamow- Teller resonance. 
Since we don't have enough data in the Kr-Tc region 
of interest, we adopt the same value as that derived in 
Ref. @ for double-beta decay emitter 76 Ge (In fact, the 
/3-decay depends on the low-lying strength distributions 
and the choice of g p h does not affect the final results too 
much) . 

The calculated half-lives are sensitive to the renor- 
malized particle-particle strength g pp . But we must 
also consider the possibility of quenching of the axial- 
vector coupling constant due to short-range correlations 
and to multi-phonon effects which are excluded out in 
QRPA calculations. Due to the lack of experimental 
data on log/t values of single decay branches in our 
mass region of interest, we take the one used in Ref. 
[Tol ] that was derived from experiment (l6[ for 150 Nd, 
gA = 0.75g.4o = 0.95, where gAo = 1-26 is the bare value 
in the vacuum. 

The relation between the calculated half-lives and val- 
ues of g pp is illustrated in Fig. [5J When g pp is increased, 
we obtain an enhanced GT strength to low-lying states, 
and therefore smaller calculated half-lives. From Fig. [21 
we find that without quenching, the half-lives are under- 
estimated, and the fitted g pp values are around zero. If 
the quenching is included, realistic values from 0.6 — 0.9 
which reproduce the half-lives of the isotopes are ob- 
tained, which agree well with the fitted g pp values of /3(3- 
decay half- lives in Ref. [T3] ■ Due to the large uncertainty 
in half-life of 100 Kr we exclude this isotope from Fig. 1. 
Another isotope which is not included in the figure is 
114 Mo, because it requires a larger model space due to 
its neutron number 72 (one more major shell should be 
added in the calculation compared with other isotopes), 
hence a much longer time is needed for calculation of the 
whole range of g pp . However, as we shall see later, the 
results for 114 Mo agree well with those obtained with the 
3pp=0.75 value we choose from the fitting. 

With the uncertainties of the choice of g pp from 0.6 — 
0.9, the errors of the half-lives vary by a factor of two in 
general. The optimal choice is g pp =0.75 from the trends 
in the upper panel of Fig. [21 In Fig. [31 we show the 
ratio between calculated and measured half-lives with 
g pp = 0.75. Following the definition from (fT2"l) . we ob- 
tain a total error of 2.1 compared to 2.3 in Ref. [l[ for 
all nuclei from Kr to Tc. If we compare our results with 
those obtained in Ref. l|] (upper panels in Fig. [3]), we 
find that we have a better agreement for even-even nu- 
clei (a total error of 1.32 compared to 2.16 for 9 nuclei 
with reasonably small error bars) not simply because we 
have adjusted the parameters for this region, but mainly 
due to the adaption of excitation energies relative to the 
ground states of the final odd-odd nuclei. This gives a 
more accurate phase-space factors which effectively re- 
duce the half-lives, and gives improved agreement with 
the experiment. While for Ref. [l|, even without the ef- 
fect of quenching, there is an over-estimation for almost 
all even-even nuclei due to the excitation energies they 



5 



10 



Cl 
X 
CD 



o 



0.1 



CL 

X 
CD 



0.1 



Even-Even 



■ Even-Odd 

▲ Odd-Event 
O Odd-Odd 



FRDM+QRPA(sep.) 



Total Error =2.16 for 9 nuclei 
-m-H 1 1 — I 



Total Error =2.19 for 22 nuclei 
' ■ ■ I " " — I 



FRDM+QR 3 A(rea.) 



Total Error =1.32 for 9 nuclei 



? o 



Total Error =2.16 for 22 nuclei 



0.01 



0.1 



0.01 



0.1 



FIG. 3: A comparison of the calculated half-lives from Ref. [l[ (upper panels) and this work (lower panels) with the experimental 
ones from RIKEN @]. The experimental errors are taken into account for several isotopes with large error bars. The total errors 
are defined in <|12[) . The left panels are comparison for even-even isotopes, and the right panels are for other types denoted by 
different symbols for the three cases of (Z, N). The three horizontal lines in each panel correspond to the ratios of 2, 1 and 0.5 
respectively. Here, sep. is the abbreviation for separable force and rea. for realistic force. 



choose. 

As an example, we show in Fig. @] the low-lying GT 
strength distributions for 110 Zr, in order to give a gen- 
eral idea of how the low-lying strength is distributed, 
and how they contribute to the decay width in deformed 
nuclei. Due to the deformations, the contributions are 
split for K — and K = ±1 parts, not only their ener- 
gies but also the strengths, this makes the distribution 
spread out. In Fig. [4j we show the position for Q value 
and the neutron separation energies. The strengths with 
the lowest excitation energies are most important to the 
decay width because of their larger phase spaces. Thus, 
nearly comparable amounts of GT strengths are located 
in the interval of S n -S2n and S^n-Q, but the /3-delayed 
two neutron emission probability P 2n is small compared 
with P n . In this sense, one needs both accurate pre- 



dictions for the strengths and their positions, and the 
advantage of realistic force is that it provides a better 
determination of the excitation energies. The effect of 
increasing g pp is that it enhances the low-lying strength 
and shifts down the excitation energies, hence reduces the 
half-lives. From the definition of P n in Refs. 0, Q, more 
low-lying strength below the neutron separation energies 
gives much smaller P n values, and vice versa. Thus, a 
comparison with the experiments for both the half-lives 
and the P n values is a good measure of how good the 
nuclear structure descriptions are. 

The agreement between experiment and theory for 
even-odd and odd-even isotopes is about a factor of two 
worse than that for even-even isotopes. The lack of 
particle-vibration coupling in the odd-mass systems at 
this region does not seem to have too much affect on 



6 




E (MeV) 



FIG. 4: As an example, the low-lying GT strength distri- 
bution (upper panel) and GT running sum (lower panel) 
for 110 Zr are shown. The three vertical lines correspond to 
S n = 4.72 MeV, S 2n = 7.79 MeV and Q = 9.79 MeV respec- 
tively, which are taken from Ref. ff|. Here, we use g pp — 0.75, 
which gives a half-life t\/2 = 37.9 ms and /3-delayed one- 
neutron emission probability P n = 10.5%. 



the final half-lives, however. This is consistent with the 
calculations in Ref. where a weak-coupling approxi- 
mation was assumed. 

The agreement between experiment and theory is 
worse for the odd-odd isotopes, and there exists a sys- 
tematic overestimation for the half-lives. This is due to 
a shortcoming of the method we use, the lack of consid- 
eration of the collectivity for even-even daughter nuclei. 
This over-predicts the energies of the final states, and 
hence the calculated phase factors are smaller than ex- 
pected. However, in spite of the shortage of the methods, 
we can still keep the error within an order of magnitude 
and for most isotopes approximately a factor of five. We 
note from [17| that the r-process path does not depend 
strongly on the beta decay properties of the odd-odd nu- 
clei due to their larger S n values. Thus, instead of im- 
proving the models for odd-odd nuclei, one alternative 
way is to simply use the average of the results for the 
neighboring odd-even and even-odd nuclei for their val- 
ues in an r-process database. 

With the above comparisons and discussions, we ex- 
tend our calculations to all of the deformed Kr-Tc iso- 
topes in the region N = 50 — 82. We make comparisons 
with experimental measurement (if available) and previ- 
ous theoretical predictions from Ref. pj . The results are 
shown in Fig. O The same set of Q- values taken from 
FRDM model as used in Ref. [l[ is adopted for the sake 
of comparison. One of the differences between our results 
and that of Ref. [l| is that the latter have added an ex- 
tra strength spreading for each of the final states. In our 
calculations of even-even nuclei the strength is already 



spread by the deformation effects, and there is not as 
much motivation for adding more by hand. But for our 
calculations of odd-mass and odd-odd nuclei where the 
collective behavior has been excluded, the transitions are 
just between the single-particle states with little spread- 
ing compared to that in Fig. [?] for even-even nuclei. In 
this case, we might get better results by adding some 
spreading in the strength. 

For most even-even isotopes, shorter half-lives by up 
to a factor of two are predicted in our calculation com- 
pared to Ref. [l[ due to the lower excitation energies for 
the final states. This behavior applies also to some odd- 
mass nuclei. For even-Z isotopes, the half-lives decrease 
with the increase of neutrons with some small stagger- 
ing behavior, but overall agreement with experiment is 
obtained. For odd-Z, there are systematic over estima- 
tions of the half-lives for the odd-odd isotopes with the 
reasons stated in Sec. II. In Ref. due to the addi- 
tional strength spreading put in by hand, they obtained 
a better agreements for these isotopes. As discussed in 
Sec. II, the r-process results are not sensitive to the half- 
lives of the odd-odd nuclei, and from a practical point 
of view it would be adequate to simply use the average 
of the calculated half-lives of the neighboring even nuclei 
for these odd-odd nuclei. 

Another observable from the experiments for some iso- 
topes is the P n value, which with accurate separation 
energies gives a measure of the Gamow- Teller strength 
distribution as shown in Fig. 01 From Fig. [51 we find 
a good agreement again for even-even isotopes from lim- 
ited data, proving that the reliability of our descriptions 
for deformed even-even isotopes in this region. In Fig. 
[5J one finds a staggering behaviours in the realistic re- 
sults for the even and odd N number neutrons especially 
for the odd-Z isotopes. The reason can be traced back 
to the treatment of the odd nuclei with the lack of the 
collectivity. This shifts the excitation energies up and 
the strength distributions are shifted systematically to 
higher energies. This behavior is more obvious for odd- 
odd nuclei where nearly all the strengths are shifted up. 
In applications of our calculations to the r-process it is 
preferable to replace the calculated P n results for odd N 
values with the average of the neighboring even N values. 



IV. IMPLICATIONS FOR THE R-PROCESS 

We have investigated how the r-process element abun- 
dances are affected by the /3-decay half-lives of var- 
ious nuclei by changing the half-lives of Moeller's 
predictions ll| for all even-even or odd-odd nuclei by one 
order of magnitude. These preliminary results agrees 
with recent r-process simulations [17j : for even-even nu- 
clei, such changes of lifetimes have tremendous effect on 
the peak formations, totally change the patterns of the 
abundance distributions. In case of odd-odd nuclei, one 
order of magnitude change in all of the half-lives results 
in essentially the same the r-process abundance pattern 



7 




58 60 62 64 66 68 70 72 74 76 60 62 64 66 68 70 72 74 76 



N 



N 



FIG. 5: A comparison among the calculated half-lives from Ref. [l|] (red), this work (blue) and measured ones (if available) 
with error bars from RIKEN Q| for Kr to Tc isotopes. The total errors are defined in Ref. Here, as before, sep. is the 
abbreviation for separable force and rea. for realistic forces. 



except for the A=150-200 mass region where the odd- 
even oscillation for elemental abundances are relatively 
changed by about a factor of two. 



We can conclude from this simple simulation that cur- 
rent accuracy of deformed QRPA method can meet the 
needs of nuclear inputs for the r-process simulation. Our 
next step is to extend the present calculations to other 
deformed regions, for example, the heavily deformed rare- 
earth elements region, where the beta-decay data is lim- 
ited, and where it is still not understood how the the 
peak of rare-earth elements is formed. The final goal is 
to calculate the /3-decay properties over the whole nuclear 
chart. It is also important to have reliable calculations 
for spherical nuclei, especially those around N — 82 that 
are important for the abundance peak around A = 130. 



V. CONCLUSION 

In this work, we investigated the /3-decay properties of 
the Kr-Tc isotopes recently measured at RIKEN. With 
the pn-QRPA taking into account of realistic forces, a 
good agreement has been obtained between the theory 
and the experiments especially for even-even nuclei, with 
an accuracy within a factor of two for most of them. 
The current calculations provide improved results for the 
beta-decay half-lives of even-even nuclei. We plan to ap- 
ply the present method to the rare-earth region of de- 
formed nuclei. We also plan to use the realistic inter- 
actions for QRPA calculations of spherical nuclei. This 
will eventually provide an improved set of predictions for 
the half-lives and P n values that can be used in r-process 
network calculations for the element abundances. 

This work was supported by the US NSF [PHY- 
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FIG. 6: A comparison among the calculated /3-delayed neutron emission probability P n values from Ref. 4] (red), this work 
(blue) and measured ones (if available) with error bars from RIKEN [j| for Kr to Tc isotopes. Here sep. is the abbreviation for 
separable force and rea. for realistic forces. The empty squares are the /3-delayed one neutron emission probability P n , while 
the solid ones are the total probability up to 3 neutron emission. 
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